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Abstract 

We found that in transition arrays of complex atomic spectra, the 
strengths of electric-dipolar lines obey Benford's law, which means 
that their significant digits follow a logarithmic distribution favoring 
the smallest values. This indicates that atomic processes result from 
the superposition of uncorrelated probability laws and that the oc- 
curence of digits reflects the constraints induced by the selection rules. 
Furthermore, Benford'law can be a useful test of theoretical spectro- 
scopic models. Its applicability to the statistics of electric-dipolar lines 
can be understood in the framework of Random Matrix Theory and is 
consistent with the Porter-Thomas law. 



1 Introduction 

Newcomb p] noticed that the beginning pages of logarithm books were more 
used than the last pages. This observation led to the conjecture that the sig- 
nificant digits of many sets of naturally occuring data describing the physical 
world are not equi-probably distributed, but in a way that favors smaller 
significant digits. For instance, the first significant digit (i.e. the first digit 
which is non zero) will be 5 more frequently than 6 and the first three sig- 
nificant digits will be 458 more often than 462. Moreover, this result is scale 
invariant. Benford provided a probability distribution function for signifi- 
cant digits. The probability that the N first significant digits dj (i = 1,N) 
are equal respectively to ki (fcj=l,9) is given by [2]: 

P{d, = h) = log 10 (l + di x 10*-*] ) . (1) 

For instance, since some accounting data satisfy this rule, it is used to 
detect fraud [3]. River lengths, mountain heigths, populations of cities [I], 
radioactive decay half lives [5], physical constants and some mathematical 
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series satisfy Benford's law. However, if many data roughly satisfy the law, 
perfect agreement is not the general case and the law is not universal. In the 
present work, focus will be put on the first significant digit, which probability 
of being equal to k is given by 

P(di =k) = log 10 (l + k- 1 ). (2) 

2 Digits of line strengths 

To our knowledge, although a great attention has been paid to the modeling 
of electric-dipolar (El) transition lines in a complex atomic spectrum typical 
of a hot plasma, it has never been checked whether the statistics of the lines 
obeys Benford's law. Most of the studies focused on the distribution of the 
values of the line strengths (see for instance [6] or [7] ) , but not on the signif- 
icant digits. Figure Q] represents the transition array 3d 6 — ► 3d 5 4p of Ge IX 
calculated with Cowan's ab initio atomic-structure code [H], which provides 
the Slater direct and exchange integrals, and the one-particle dipole-moment 
integrals evaluated from Hartree-Fock wavefunctions. Figure [2] displays the 
number of lines of that transition array for each value of the first significant 
digit calculated from our approach and predicted from Benford's law. The 
total number of lines is 3245. It is clear that the two distribution of lines 
are very similar. The same property has been observed on other transition 
arrays for different elements. Table [T] displays the fraction of lines per sig- 
nificant digit for four transitions: 3d 6 — > 3d 5 4p (3245 lines) for Ge IX (Ti), 
3d 4 -> 3d 3 4f (2825 lines) for Ta LII (T 2 ), 3d 2 4s4p -> 3d4s4p 2 (2722 lines) 
for Br XIV (T 3 ) and 3d 2 4s4p 3 -» 3d4s4p 4 (8231 lines) for Br XII (T 4 ). It 
is worth mentioning that even for transition arrays with a much smaller 
number of lines, the law still applies very well. 

3 Testing spectroscopic models 

In a hot plasma, the total spectrum consists of a huge number of transition 
arrays. Therefore, the total number of lines can be immense. Usually, when 
the lines coalesce due to the physical broadening processes (natural width, 
electron impact, Stark effect and Doppler effect), the distribution of lines 
in intermediate coupling is modeled by a continuous (Gaussian) distribu- 
tion, named UTA (Unresolved Transition Array) which average energy and 
variance are evaluated following the work of Bauche et al. [9|. However, 
in that case some information is lost, and random calculations can help 
"re-constructing" approximately the resolved distribution of lines ensuring 
known physical rules [TO] . 

The RTA (Resolved Transition Arrays) approach developed by Bauche 
et al. [H] and recently improved by Gilleron et al. [12] enables one to 
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Table 1: Fraction of lines which first digit is k for transitions Ti, T2, T3 and 
T4 compared to the values predicted by Benford's law. 



simulate El lines of a transition array avoiding direct diagonalization of the 
atomic Hamiltonian. The strength of a line and the energies of the lower and 
upper levels are picked up at random satisfying some specific constraints: 
preservation of the number of lines, of the total strength and weighted vari- 
ance of the transition array, as well as of the unweighted variances of levels 
of lower and upper configurations. In that approach, the distribution func- 
tion also includes a correlation between the energies and amplitudes of the 
lines (propensity law). Therefore, Benford's significant-digit law can be a 
performing test of the validity of the correlation rules included in such ap- 
proaches. 

4 Interpretation 

The fact that El lines satisfy Benford's law indicates that the distribution 
of digits reflects the symmetry due to the selection rules. If transitions 
were governed by uncorrelated random processes, each digit would be equi- 
probable. 

Benford's law is still not fully understood mathematically. However, 
some key elements are known since the middle of the nineties. For instance, 
Benford's law holds for any rescaling of the data, since it does not depend on 
any particular choice of units [21J. It is the only scale-invariant law referring 
to digits since the logarithmic distribution is the unique continuous base- 
invariant distribution [22j . More precisely, if probability distributions are 
selected at random and random samples are taken from each of these distri- 
butions so that the overall process is base invariant, then the significant-digit 
frequencies of the sample will converge to the logarithmic distribution. In 
other words, the more diversified the probability distributions are, the bet- 
ter the data sets fit Benford's law. The understanding of why systems with 



3 



many interacting particles spontaneously organize into scale-invariant states 
is a difficult task of statistical physics. 

Moreover, Benford's law applies if the system is governed by multiplica- 
tive processes [23]. Indeed, a random multiplicative process corresponds to 
an additive process (dynamical description of a Brownian process) in a loga- 
rithmic space. In Wigner's Random Matrix Theory |14|. I15j . the Hamiltonian 
is defined in the Gaussian Orthogonal Ensemble (GOE) by an ensemble of 
real symmetric matrices which probability distribution is a product of the 
distributions for the individual matrix elements H^, which are considered 
as stochastic variables. The variance of the distribution for the diagonal 
elements is twice the one for the off-diagonal elements. The line strength S 
is proportional to | < i\D\j > | 2 , where D is the dipolar operator and \i > 
and \ j > are eigenvectors of the Hamiltonian. Therefore, the line strengths 
involve quantities such as products of H^i and one has 

f = C, (3) 

where S and S' are the strengths of two lines belonging to the transition 
array and £ is a stochastic variable. Eq. ([3]) can be written 

log S = log C + log S' (4) 

The central-limit theorem states that the probability distribution that 
the value of the n th strength is S will be Gaussian with a variance cx n 1 / 2 . 
In the infinite limit, the distribution will approach the uniform one, charac- 
terized by constant K. Therefore, one has 

f P{logS)d(\ogS)=K (5) 

The probability P that the first significant digit d\ of S is k in base 10 
is given by 

f k+1 dS I f 10 dS , / 1\ 

which is exactly Eq. ([2]). This heuristic argument relies on the idea that 
fluctuations are governed by multiplicative processes involving a stochastic 
variable. The matrix elements of the Hamiltonian are correlated stochastic 
variables and the product of such variables leads to Benford's logarithmic 
distribution of digits. 

Another argument can be invoked to understand the mathematical fun- 
dation of Benford's law. Porter and Thomas have shown that the amplitudes 
of the lines between all the levels of two random matrices obey a Gaussian 
distribution [TBI HZ]- ^ ne strength being the square of the amplitude, its 
distribution is [6j [18] 
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where L and < S > are respectively the number of lines and the average 
value of the line strength in a (J, J') set. When numbers are taken from an 
exponential distribution, they automatically obey Benford's law. Therefore, 
when the exponential term dominates in D(S), which is often the case except 
close to the origin (i.e. for very weak lines), Benford's law applies. 

The Random Matrix Theory contains approximate symmetries, which 
are not sufficient to describe the vicinity of Russell-Saunders and jj cou- 
plings |19| 120] . In the model of Wilson et al. |13j . diagonal terms are 
calculated in a pure coupling using Cowan's code mentioned above and off- 
diagonal elements are populated statistically beyond the GOE according to 
a bi-Gaussian distribution function where elements are correlated. They 
observed a disproportionally large number of off-diagonal elements of small 
amplitute. However, we found that even when the Random Matrix The- 
ory is expected to be inappropriate, i.e. close to a pure coupling, the line 
strengths still fit Benford's law. This is due to the fact that even if the 
number of weak emerging lines is important (when the term -7= dominates 
in the Porter-Thomas law ©), the number of decades is sufficiently large so 
that the weight of those weak lines is negligible in the statistical occurence 
of digits. 

5 Conclusions 

The results of this study are following. First, the distribution of lines in 
a given transition array follows very well Benford's logarithmic law of sig- 
nificant digits, which was, to our knowledge, never observed and rather 
unexpected. 

Second, this indicates that the correlations due to the selection rules 
manifest themselves in the distribution of digits, even in intermediate cou- 
pling, and that such a symmetry results from a superposition of different 
probability laws. 

Third, computational methods avoiding diagonalization such as the Re- 
solved Transition Arrays model, which ensures known physical rules (propen- 
sity law, correlation laws), can be improved and/or tested by imposing the 
constraint that the statistics of the lines must verify Benford's law with a 
high accuracy. 

Finally, since Benford's law can be explained in terms of a dynamics gov- 
erned by multiplicative stochastic processes, the Random Matrix Theory is 
probably an interesting pathway for the calculation of large atomic-dipole 
transition arrays and Benford's law can help clarifying the existence of dif- 
ferent classes of stochastic Gaussian variables. Moreover, Benford's law is a 
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signature of the Porter-Thomas distribution of lines. 
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Figure 1: Array of lines for the transition 3<i 6 — > 3(f 5 4p of Ge IX. 
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Figure 2: Fraction of lines versus first significant digit for the transition 
3d 6 — ► 3d 5 4p of Ge IX. The total number of lines is 3245. 
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